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Abstract. - The conductance in two-dimensional (2D) normal-superconducting (NS) systems is 
analyzed in the limit of strong magnetic fields when the transport is mediated by the electron- 
hole states bound to the sample edges and NS interface, i.e., in the Integer Quantum Hall Effect 
regime. The Andreev-type process of the conversion of the quasiparticle current into the superflow 
is shown to be strongly affected by the mixing of the edge states localized at the NS and insulating 
boundaries. The magnetoconductance in 2D NS structures is calculated for both quadratic and 
Dirac-like normal state spectra. Assuming a random scattering of the edge modes we analyze 
both the average value and fluctuations of conductance for an arbitrary number of conducting 
channels. 
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Introduction. — Andreev transport phenomena, i.e., 
transport effects associated with the conversion of elec- 
tjrons into holes, are known to determine the distinctive 
features of a wide class of hybrid structures consisting of 
tihe normal (N) and superconducting (S) metal parts (see, 
e.g., [T] and references therein). Applying rather high 
ijiagnetic fields one can drastically affect the physics of 
these Andreev-type effects due to a strong modification 
qf the transport mode structure [5]-[S] which is typical for 
the systems in the Integer Quantum Hall Effect (IQHE) 
r'egime. Provided the radii of the cyclotron orbits in the 
normal part of the system become less than the mean free 
path the transport appears to be determined by the waves 
bound to the sample edges. Depending on the momenta of 
these waves and quasiparticle charge the transport modes 
are localized near the different edges and, thus, the wave- 
functions of the incoming and outgoing particles appear to 
be spatially separated. Thus, the magnetic field destroys 
the basic backscattering property of the standard Andreev 
reflection. 

It is the goal of the present work to suggest a general 
theoretical description of the Andreev transport mediated 
by the edge states in the IQHE regime. We consider here 



an exemplary two-dimensional (2D) NS system shown in 
Fig. [1] Such type of mesoscopic junctions based on a 2D 
electronic gas (2DEG) or gapless 2D semiconductors like 
graphene are in the focus of current experimental and the- 
oretical research PHTU]. To elucidate our main results we 
start here from a qualitative description of the transport 
mediated by the edge states. An electron injected from 
the normal conductor goes to the superconductor through 
an edge state "a". At the "ab" -corner it transforms into 
two types of hybridized electron-hole states at the bound- 
ary "b" with the probabilities ti and (1 — ti), respectively 
(see the inset of Fig. [1]). Similarly to the situation at 
the "ab" -corner each of electron- hole quasiparticles trans- 
forms with the probabilities T2 and (1 — T2) into an elec- 
tron and a hole which return to the normal lead through 
the edge states "c". Without mode mixing (ri_2 = or 
1) each initial electron (hole) state at the boundary "a" 
completely transforms into the final electron (hole) state 
at the boundary "c" and the probability of the electron- 
hole conversion is zero, so that the total conductance G 
vanishes. Taking account of the quasiparticle mode mixing 
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at the corners we find: 

G = 2GoiV[ri(l-T2)+T2(l-Ti)] , (1) 

where Gq = e^/Trfi is the conductance quantum, and N is 
the number of propagating electron edge states. Thus, it 
is the mode mixing which is responsible for the conversion 
of the quasiparticle current into the supercurrent outgoing 
from the NS boundary. 

Generally, the solution of the problem of coupling be- 
tween the edge states near the corner is rather complicated 
and depends on the details of the system geometry. There- 
fore only limiting cases have been previously considered: 
(i) quasiclassical limit with the large number of the edge 
states, when the problem can be treated considering the 
particles and holes on the cyclotron orbits skipping along 
the surface [3]; (ii) quantum limit when the number of 
the edge states is of the order of unity [3J[S]. The first 
case corresponds to rather large Fermi energies comparing 
to the Landau level spacing. The corresponding oscillat- 
ing behavior of the conductance of the 2D NS junction 
vs magnetic field, junction width and/or Fermi level has 
been analyzed in detail in Rcfs. 

In this Letter we analyze the magnetoconductance be- 
havior in the NS structures for an arbitrary number of 
transport modes taking into account the mixing of the 
edge states near the corners. Adopting a phenomenologi- 
cal description of the mode mixing problem based on the 
transfer matrix approach we find a simple expression for 
the conductance in the quantum limit. Such model sug- 
gests a simple explanation of the oscillatory phenomena 
mentioned above and brings out their dependence on the 
mode coupling parameters. To find the conductance for 
an arbitrary number of quantum channels we assume the 
mode mixing at the corners to be random in the sense 
that the appropriate scattering matrices arc uniformly dis- 
tributed (see below for details). Such approach allows us 
to find universal expressions for the average conductance 
and its fluctuations which depend only on fundamental 
constants and number of transport channels: 

(G) = GoA , (2) 

-c^aA^^^^Go^II^. (3) 

In the equation ([2]) contribution of each electron trans- 
port mode to the conductance equals to the conductance 
quantum Go instead of 2Go in ballistic NS-junctions with- 
out magnetic field. Such conductance reduction is caused 
by the levelling of outgoing (along the "c"-edge) electron 
and hole probabilities in the limit of strong disorder, which 
saturate at the value 1/2. Such levelling is analogous to 
the one observed in numerical simulations in [5] for the 
tight-binding model with the disorder in on-site energies. 

Basic equations and edge state spectra. The 

spectra of quasiparticle edge states can be found using 



the Bogolubov-de Gennes (BdG) equations written for 
electron-like (u) and hole-like (w) parts of the wave func- 
tion ^ ~ {u, i)): 

Here T is the time-reversal operator, A is the gap oper- 
ator, and the energy E is measured relative to the Fermi 
level /i. Note that we neglect here the Zeeman shift of the 
quasiparticle spectra caused by the interaction of magnetic 
field with the true electron spin. 

For 2DEG the single particle hamiltonian H takes the 
Schrodinger form H{A) = (p — |A)^/2m, where p = 
—ih{d/dx;d/dy) is the momentum operator in the x-y 
plane of the 2D system, and A is the vector potential cor- 
responding to the magnetic field B perpendicular to the 
system plane. In graphene there appear 2 sublattice (pseu- 
dospin) and 2 valley (isospin) degrees of freedom and in the 
"vallcy-isotropic" basis the hamiltonian could be written 
as follows: H = vpTo (g) cr (p — |A) (see, e.g., [1]). Here 
wf is the Fermi velocity, ct^ and are the Pauli matri- 
ces acting in the sublattice and valley spaces, respectively, 
and (To, To are the 2x2 unit matrices. 

The boundary conditions and corresponding edge state 
spectra at the boundaries with isolator and superconduc- 
tor have been previously studied for 2DEG |2] as well as 
for graphene [TOTHT^ . The boundary condition at the NS 
interface with 2DEG couples electron u and hole v parts 
of wave function and in quasiclassical limit takes the usual 
form V ~ e^^^u, where (3 = arccos(e/A). The wave func- 
tion at the 2DEG insulating edge vanishes 'I' = 0. Accord- 
ing to Akhmerov and Beenakker [J] the graphene-isolator 
boundary conditions 

* = (i^, f) ® (n_L,CT)* (5) 

are crucially determined by the isospin vector v while the 
resulting quasiparticle spectrum depends on the vector 
nj_. Here the unit vector n_L should have a zero projec- 
tion on the direction normal to the graphene edge. The 
graphene-superconductor (GS) interface boundary con- 
dition doesn't depend on valley degree of freedom and 
for subgap energies it could be written as follows: v — 
^-tfi(n-cr)^^ Thus, the edge state spectrum is valley degen- 
erate. 

Taking the case of a homogeneous boundary and choos- 
ing the gauge with the vector potential parallel to the 
boundary one can find a set of spectral branches En{k\\) 
vs the conserved momentum component fc|| along the in- 
terface. Here we introduce an integer index n enumerating 
the branches. Thus, each insulating edge ("a" or "c") sup- 
ports 2 A propagating edge modes: A electron- like modes 
and A hole-like modes. The NS interface also supports 2A 
propagating (valley degenerated in graphene case) modes 
with mixed electron- hole wave functions. 

Mixing of the edge modes. Transfer matrix ap- 
proach. — Considering the transport mediated by the 
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edge states we use a simple phenomenological model based 
on the transfer matrix approach. We introduce a trans- 
fer matrix S which couples the quasiparticle edge waves 
propagating along the "a" and "c" boundaries. This ma- 
trix calculated for the states at the Fermi level is known 
to determine the linear transport characteristics at zero 
temperature [I1[I1|- It is important to note here that at 
the Fermi level the matrix S appears to be simultaneously 
a scattering matrix coupling the incoming and outgoing 
electron-hole waves. Indeed, in this case all the states 
propagating along the interfaces in Fig. [1] have the same 
sign of the group velocity dEn/dk\\ [2111] and, thus, all 
quasiparticle fluxes are flowing clockwise for a chosen mag- 
netic field direction. 

The BdG equations dH) are known to be invariant with 
respect to the transformation converting electrons into 
holes and changing the sign of energy E and, as a conse- 
quence, all the edge modes can be divided into two groups 
connected by this transformation. For the NS interfaces 
we denote these groups as (7+ and g- while for the bound- 
aries with an insulator ("a" and "c") these groups just 
coincide with pure electron [ua and Uc) and hole {va and 
Vc) waves. Each of hybrid states g± propagating along 
the ballistic NS boundary "b" of the length L acquires a 
phase factor e^'^*'"'"^, where kn is the momentum satisfy- 
ing the equation En{kn) = for the spectral branches at 
the NS boundary. The scattering processes at the corners 
"ab" and "be" (see Fig. [T]) could be described by uni- 
tary transfer matrices Ti 2, which couple the incident and 
transmitted quasiparticle waves: 



.9+ 
.9- 



Ti 



9+ 
9- 



(6) 



where Ua, Va, Uc, Vc, g+, g- arc the sets of wave ampli- 
tudes corresponding to the solutions of BdG equations at 
the Fermi level. The unitarity of the matrices Ti^2 is a con- 
sequence of the quasiparticle current conservation. The 
scattering matrices Ti^2 can be conveniently presented in 
the four-block form: 



t+e t+h 
t-e t-h 



To = 



th+ 



4- 



(7) 



The total transfer matrix S can be written as a product of 
three matrices describing subsequent scattering and prop- 
agation processes discussed above: 



S 



She 



Seh 
Shh 



T2ALT1 , 



(8) 



,6"''="^) is a 



where A = diag (e**^!^ 
diagonal transfer matrix of phase factors acquired at the 
NS boundary. The blocks See, Se/i, She, Shh of the matrix S 
describe the scattering between the edge states at the "a" 
and "c" boundaries. The zero-temperature conductance 
is given by following expression [14] : 



G = GoSp /-s 



+ stiJeh) - 2GoSp (s+ Seh) . (9) 



Magnetoconductance in NS junctions. Quantum 

limit. — At sufficiently small Fermi energy /i each edge 
supports only two propagating modes and each block of 
scattering matrices becomes a single complex amplitude. 
It is convenient to parametrize these amplitudes as follows 



t+e = ^/^^e,'(^^+''^/2) , t_e 



Tie 



j(tpi-ei/2) 



te+ = ^/T^e^(^^+^^/2) , ie- - V^t 



.1(^2-62/2) 



(10) 
(11) 



=^ ±t 



±e. 



th^ 



±t*^, where ti is the probability 



that electron at the boundary "a" scatters into the second 
type of hybrid modes, T2 is the probability that the first 
type hybrid mode scatters into the hole at the boundary 
"c" . The matrix of phase factors takes the form: A = 
diag (e**''^'^, e"**"'^^) , where fci is the momentum value at 
which the spectral branch at the NS boundary crosses the 
Fermi level. Omitting the calculation details we present 
here the final expression for the 2DEG - superconductor 
junction conductance: 



G = 2Gq{p2i +P12 - cos 7) 



(12) 



where p„„ = r„(l - r^), 7 = 2kiL + 9i + 92- In the 
symmetric geometry of Fig. [TJ i.e., when the "ab" and 
"be" corners have the same shape, one can expect the 
appearance of an additional symmetry of scattering ma- 
trices describing the mode mixing: T2 = 2f (pi2 = P2i)- 
The expression for conductance in this case can be further 
simplified: 



G = 4Gor(l-r)(l-cos7) 



(13) 



where ti = r2 = r. The conductance reveals an oscil- 
lating behavior vs the junction width L which is, in fact, 
a consequence of quantum mechanical interference of the 
edge waves propagating along the NS boundary. The ex- 
pression ()13p is in good agreement with the qualitative 
arguments in the introduction: the absence of mode mix- 
ing corresponding to the limits r = or r = 1 causes a 
complete suppression of the charge transport. 

Considering a quantum limit for GS junctions we need 
to emphasize two important distinctive features. First, 
the momentum of the zero energy mode (i?i(fci) = 0) at 
the GS boundary appears to vanish [ki = 0) and both 
states [g+ and g-) become degenerate. Second, the scat- 
tering matrices Ti 2 crucially depend on the isospin degree 
of freedom. Following Ref. |4] we introduce the isospin op- 
erator eigenvectors' basis {vm,T) \±Vm) = ± jit'm), where 
vi,2 are the isospin vectors characterizing the "a" and "c" 
boundaries (see the boundary condition ([S])). The scatter- 
ing matrices take the form: 

fi = {vi\t+, + \-vi) {Vl\t^,,+ 

+ {-vi\t+h + {-yi\t-h (14) 

T2 = |l^2) (1^2! te+ + \-iy2) {l^2\th+ + 

+ W2){-l^2\te- + 1-1^2) {->^2\th- (15) 
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Introducing the notation x — {^1-1^2) for the angle be- 
tween the isospin vectors one can get the conductance of 
the system in the form: 

G = 2Go [sin^ | + (p2i + P12) cosx 

^ 2v'Pi2P2i(cos 61 cos 02 — sin 61 sin 62 cos x) 

- sinx((l - 2Ti)y/p^sin0i +(1 - 2r2)^/pITsin 6*2)] . 

(16) 

Assuming a symmetric geometry of Fig. [T]we put T2 = 
and find: 



G = 2Go (2 VP • cos f sin 61 - (1 - 2t) sin f ) 



(17) 



where 9^ = 62 = 0, p = r(l - r). Contrary to the 2DEG 
case the conductance does not depend on the junction 
width L which is a natural consequence of zero phase ac- 
quired by the waves propagating along the GS boundary 
in the two-mode limit. Another new feature specific for 
the case of graphene is that the additional mode mixing 
occurs due to the mismatch of the isospin directions at 
different insulating boundaries. Note that neglecting the 
intervalley scattering at the corners we put r = or r = 1 
and get the limit considered in Ref. [1]. 

Magnetoconductance in NS junctions. Random 
matrix theory. — Considering the charge transport me- 
diated by a large number of edge states it is natural to 
expect that the conductance will be given by the sum of 
phase factors: 



i{ki — kj)L 



(18) 



where the hermitian matrix is determined by the trans- 
fer matrix parameters. The oscillating behavior of conduc- 
tance vs L becomes more complicated than in a two-mode 
limit and is generally characterized by a set of incommen- 
surable periods. Previously these oscillations have been 
predicted on the basis of quasiclassical method for quasi- 
particlcs moving along skipping cyclotron orbits [3]. In 
real experimental situation the mode interference and the 
corresponding oscillations should be, of course, smeared 
due to the effect of sample imperfections, e.g., roughness, 
etc. Here we suggest a phenomenological approach to treat 
the problem taking account of these effects for an arbitrary 
number of modes. We model the scattering caused by the 
sample imperfections introducing random transfer matri- 
ces Ti^2 and applying standard methods for calculation of 
the ensemble averages [1]. 

The symmetry properties of the transfer matrix Tm al- 
low us to introduce a polar decomposition (cf. [T51IT5] ): 
Tm = d„e„/3„,, where 





A* 


) , An = 1 


^ Bm 

V 




- Tm V 






fr^i Vl 









and the unitary matrices Am,Bm characterize the scat- 
tering phase shifts while the diagonal matrix fm = 
diag(rm', ■ . • ,T^n'^) consists of the eigenvalues t{"^ (^'2"^) 
of the matrix i'^J.-e {it-^e-)- The values Tm'' give us 
the probabilities of transitions between the modes. Any 
scattering 2A^ x 2A^-matrix X invariant under electron- 
hole converting transformation of BdG equations belongs 
to the compact symplectic group Sp(A'^) = {X2Nx2N ■ 
XX+ ^i,JX + X^J = 0}, where 



J 



/ 

-/ 



Our further calculations are based on the simplest as- 
sumption about the distributions of the random transfer 
matrices: we consider the random unitary matrices Tm 
uniformly distributed in the compact symplectic group 
Sp(A^). 

The uniform distribution of an element A of the compact 
group is defined with respect to a measure dfi^A) which is 
invariant under multiplication: dfi{A) = dii{UAV) for ar- 
bitrary elements C/, V belonging to this group. This mea- 
sure is known as the "invariant measure" or "Haar mea- 
sure" [T] . Note that under such assumption all distinctive 
features of the graphene case associated with the isospin 
degree of freedom do not reveal in the averages and, thus, 
further expressions are valid for both types of junctions 
under consideration. 

Analogously to Ref. we derive a full distribution of 
Tm for arbitrary N and the statistics of the eigenvalues 

dfl{fm.) = P{{Tm})Y[dT^Ufl{Am)dKBm) , (19) 

a 

where P{{Tm}) ^ C Y[ hm —Tm^l is the joint probability 

distribution of the {Tm} values, dfj,{Am) [d^{Bm)] is the 
invariant (Haar's) measure on the unitary group for matrix 
Am. [Bm] and C is a normalization constant. 

The conductance averaged over the matrices Am., Bm 
takes the form: 



{G)^ B = 2GoN [n(l - f2) + f2{l - fi)] 



(20) 



where fp = Sp(fp)/A^. This expression gives us a gener- 
alization of the Eq. (|12p averaged over the phase 7 and 
written for an arbitrary number 27V of transport modes. 
After averaging over the eigenvalues Tm^ we find the ex- 
pressions and (O for the conductance and its square 
deviation. The ensemble average conductance is propor- 
tional to the number of channels and doesn't depend on 
junction width L. The square deviation increases with 
the growing channel number and saturates at the univer- 
sal number Go/2\/2. 

Changing the applied magnetic field B one could control 
the number of edge modes which decreases stepwise with 
the increasing field. 
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Thus, the dependence of conductance vs the inverse 
magnetic field B^^ reveals a series of equidistant steps 



G 



2d 



Go 



mcii 
ehB 



Go 



Go 



C/i 



2vlehB 



for 2DEG and graphene respectively. Here we introduce 
the notation [...J for the integer part. 
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Fig. 1: (Color online) The junction between the superconduc- 
tor and the chiral metal [edge states in high mobility 2D metal 
in the IQHE regime]. Solid (dotted) lines correspond to the 
electron (hole) edge states at "a" and "c" insulating bound- 
aries, and wavy lines correspond to the hybridized electron- 
hole modes at the superconducting boundary "b". Inset: Edge 
mode scattering from insulating edge "a" to superconducting 
edge "b" with probabilities ri and (1 — ti). 
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